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We have generalized the application of cumulant expansion to ferrimagnetic systems of large 
spins. We have derived the effective Hamiltonian in terms of classical variables for a quantum 
ferrimagnet of large spins. A noninteracting gas of ferrimagnetic molecules is studied systematically 
by cumulant expansion to second order of (Js/T) where J is the exchange coupling in each molecule, 
s is the smaller spin (Si,S2) and T is temperature. We have observed fairly good results in the 
convergent regime of the expansion, i.e T > Js. We then extend our approach to a system of 
interacting ferrimagnetic molecules. For one dimensional nearest neighbor interaction we have 
observed that the correlation of more than two neighboring sites is negligible at moderate and 
high temperature behavior. Thus the results of a single molecule can be applied to the chain of 
interacting molecules for temperatures greater than classical energy scale, i.e T > JS1S2. Finally 
we will discuss the effect of spin inhomogeneity on the accuracy of this method. 

PACS numbers: 75.40.-s, 75.10.Hk, 75.10.Jm 



I. INTRODUCTION 

Ferrimagnets are systems composed of different spins, 
mostly of two types (S\,S2)^- Many of ferromagnetic 
materials are essentially ferrimagnets 3 . It is a special 
class of spin models with much attraction in low dimen- 
sion where their universality class is completely differ- 
ent from homogeneous spin models^. Antiferromagnetic 
chain of spin s models are gapless (gapped) for half- 
integer (integer) spins 5 . However ferrimagnetic chains 
behave differently. The low energy spectrum of a fer- 
rimagnetic chain (1,1/2) is gapless with ferromagnetic 
order and a gapped band above it which has antifer- 
romagnetic properties^. This effect causes a crossover 
from ferromagnetic to antiferromagnetic behavior at fi- 
nite temperatures^. 

Spin wave theory gives an explanation for low temper- 
ature physics of ferrimagnetic chains 10 . However, it is 
valid for temperatures smaller than the classical energy 
scale (T < JS1S2). It is our task to obtain the physics 
of ferrimagnets systematically at moderate (where quan- 
tum corrections are important) and higher temperatures. 
Cumulant expansion (CE)ii*i2ii£ is our approach to this 
problem. Recently this approach has been implemented 
to study finite temperature behavior of homogeneous 
large spin systemsi^i^. It has been observed that cu- 
mulant expansion converges in a region T > Js which 
is wider than corresponding one in the high temperature 
series expansion and even overlaps the validity regime of 
spin wave theory for homogeneous spin systems. 

We have generalized the application of cumulant ex- 
pansion to ferrimagnetic models. In this approach we 
obtain a quasiclassical Hamiltonian for a (Si, S2) system. 
The effective Hamiltonian which is a function of classical 
variables takes into account the quantum corrections sys- 
tematically in the order of (Js/T) where s is the smaller 
value of (Si, S2). The partition function is calculated by 
using the effective Hamiltonian to get physical proper- 
ties. Apart from the effective Hamiltonian for a general 



(Si,S2) ferrimagnets we found that the n-th order ef- 
fective Hamiltonian contains correlation of (n + 2)-sites. 
Moreover, for nearest neighbor interaction in ferrimag- 
netic chain, the first order effective Hamiltonian is com- 
posed of 2-sites correlation. Consequently the results of 
CE for a single molecule of {Si, 83) is the same as one 
dimensional interacting molecules. We have also found 
that in the convergence region, the gas model is a good 
representation of a chain of interacting molecules up to 
2nd order expansion. 

It is remarkable that in the special case Si = S2, our 
results match the recent studies on homogeneous spin 
chainsiii^. 

The outline of this paper is as follows. In scctionHTIwc 
have introduced the cumulant expansion in spin coherent 
state representation. The effective Hamiltonian and par- 
tition function of two sites gas model are derived in sec. 
11111 In section HVI we have derived the effective Hamilto- 
nian of ferrimagnetic chain and then the physical proper- 
ties. Results and discusion have been demonstrainted in 
sec0 where we have presented some physical quantities 
of a chain with nearest neighbor interaction. 



II. CUMULANT EXPANSION IN SPIN 
COHERENT STATE REPRESENTATION 

The cumulant of N operators or classical variables (Aj) 
is defined asii*i£ 



(Ai...A n ) c 
d d 
dXi'" dX n 
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where (■ • ■) means averaging over a classical distribution 
function or a quantum stated. It is easy to show that in 
case of all A equal to A, one will arrive at the following 
equality, 



ln(e AA ) = (c - 1) 



(2.2) 
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To do so, multiplying by A" /n! and summing over n = 
1, 2, oo is performed. Let us show the cumulant of few 
cases in terms of usual averaging. 

(A) c = (A), (A 1 A 2 ) C = (A 1 A 2 )-(A 1 )(A 2 ) 
(A 1 A 2 A 3 ) C = (A X A 2 A 3 ) - (A 1 )(A 2 A 3 ) C - (A 2 )(A x A 3 y 
- (A 3 )(A 1 A 2 ) C -(A 1 )(A 2 )(A 3 ) (2.3) 

It is obvious from Eq. 12.3fl that the cumulant of two or 
more operator contains correlations of them. It is a char- 
acteristic property of quantum mechanics that the expec- 
tation value of a product of two operators is different from 
the product of expectation values ((AiA 2 ) — (Ai)(A 2 ) ^ 
0, whereas the equality is valid in classical limit). This is 
the correlation of two operators which is inherited from 
the quantum state in which the expectation is calculated. 
So we expect to obtain a quasiclassical description by us- 
ing cumulants. In this respect we are going to obtain an 
effective (quasiclassical) Hamiltonian in terms of cumu- 
lants which contain quantum corrections. 

The partition function (Z) of a spin system can be ex- 
pressed in the basis of spin coherent states |n); i.e. states 
with the maximum-spin projection on the axis pointing 
in the direction of unit vector n . The basis of coher- 
ent states are overcomplete, so it contains all quantum 
states. The classical state of a spin is achieved in the 
limit S — > oo. The unity operator in this representation 
is 
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Thus the effective Hamiltonian (Ti) can be considered as 
a systematic expansion in terms of cumulants of powers 
of H. The first term (ft (0) ) is the pure classical contribu- 
tion and higher orders are responsible for quantum cor- 
rections. In fact this is an expansion in powers of (Js/T) 
where s is the smaller spin of (Si,s 2 ) in a ferrimagnetic 
chain, (it will be shown in next sections). This expan- 
sion is justified whenever T > Js. To the best of our 
knowledge, a high temperature expansion for mixed spin 
systems (Si, s 2 ) is still missing in the literature. However 
to make a qualitative comparison we consider the homo- 
geneous cas o 14 ! 15 (Si = s 2 = S). In this case the conver- 
gence region of CE; i.e T > JS, is larger than the typical 
high temperature expansionist which is T > JS 2 . As 
we will see in next sections CE can also be interpreted as 
an expansion of 1/s. So for large spin ferrimagnets the 
range of convergence is wide. 

III. TWO SITES GAS MODEL 
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1 = ^T J dn\n)(n\. (2.4) 

The trace of an operator in a single spin problem can be 
written as: 

tr(A) = 2j(m|A|m) = — / dn^(m\A\n) (n\m) 

m m 

= J dn(n\A\n). (2.5) 

So the partition function of a system of TV-spins defined 
by the Hamiltonian H , is as the following, 

Z = J nl 1 (^-±)dn l (m---n N \e-^\m---n N ) 

(2.6) 

where (3 — 1/T (choosing k& = 1). If we define 

e-^ = (n 1 ...n J v|e-^|n 1 ...n JV ) (2-7) 

then the partition function is the same as a classical one. 
Using Ea. H2.2fl . the effective Hamiltonian is expressed in 
terms of cumulants. 

W=(ni--- n N \l - e-^|m ■ • • n N ) c (2.8) 

The above equation can be expanded in the following 
form by a Taylor expansion. 

m = P{H) c -^(HH) c + ^(HHH) c + --- 

= 13(H {0) +H (1) +H i2) + ■■■) (2.9) 



Let us first consider an ideal gas-noninteracting-of fer- 
rimagnetic molecules (gas model) where each molecule 
composed of two spins (Si, s 2 ). The interaction between 
spins in all molecules is either ferromagnetic or antifcr- 
romagnetic given by the following Hamiltonian. 

H = TJSi-s 2 , J>0 (3.1) 

where Si and s 2 are spin operators of size Si and s 2 
respectively. The (— ) sign stands for ferromagnetic (F) 
and (+) for antiferromagnetic (AF) interaction. There 
are two reasons to study the gas model. Firstly, it is ex- 
actly solvable, so we will compare the results of cumulant 
expansion with the exact ones to infer the accuracy of our 
method. Secondly, we would like to address the connec- 
tion of this model to a chain of interacting ferrimagnetic 
molecules (next section). 

The eigenstates of Eq. l|3.1|l are labelled by total spin 
St — \Si — s 2 \, . . . , Si + s 2 . Each state is 2S t + 1 degen- 
erate whose eigenvalue is E n = T^[St(St + 1) — Si (Si + 
1) — S2(s2 + !)]• I n the incoming calculation of cumu- 
lants we define s = min{Si, s 2 } = s 2 and the following 
parameters, 

lu= — >1. J = Js 2 . (3.2) 

s 2 

In order to achieve the cumulants of Eq. I|2.9|l it is conve- 
nient to express the spin operator on each site in the co- 
ordinate system with the z axis along the coherent state 
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vector rii = nf . 



following expression. 



S»= n TSi\ n±E( ni Tm 9 )/2 (3.3) 

OLi— Z,± 

Where and n y are appropriate transverse basis vec- 
tors. Expanding (Si -S2) in terms of spin components de- 
fined in Eq. I|3.3|) the only nonzero term in coherent state 
representation is ((Sfnx) ■ (s^^)) = (ni ■ r^^Sfsf). So 
the classical contribution to the effective Hamiltonian in 
Ea. H2.9|) is: {H) c = =F Jw(ni • n2). Quantum corrections 
to the zeroth order approximation (classical) appear in 
the remaining terms. Keeping the first two corrections 
one should calculate the following expressions. 



(HH) C = J 2 ((S 1 -s 2 )(S 1 -s 2 )) c 



(3.4) 



{HHH) C = TJ 3 ((Si •s 2 )(Si •s 2 )(Si • s 2 )) c (3.5) 

Using Ea. l|3.3|l there are 81 terms in (HH) C among which 
three are nonzero. Table (QJ shows the nonzero terms of 
four operators cumulant. A similar calculation has been 
done for (HHH) C to find the nonzero terms. In this 
case each cumulant contains the product of 6 operators 
where only 15 of them are nonzero. Summing up, we will 
arrive at the effective Hamiltonian up to second order 
corrections (0(l/s 2 )). It is written in terms of classical 
variables where quantum effects has been considered. 

H = J[T^(ni-n 2 ) 

-i[^)^ + l)(l-( ni .n 2 ) 2 ) 



1 / (3 J 



co(l - ni • n 2 ) 2 



(uj 2 - (3w 2 + uj 3 + w)(n! • n 2 ))] (3.6) 



The first term shows the classical contribution which sim- 
ply represents the energy of coupled classical spins whose 
lengths are S± and s 2 . Quantum corrections have a non- 
Heisenberg form and their structure becomes more and 
more complex with increasing order. These corrections 
are important in the intermediate temperature where 
classical fluctuations are not strong enough to suppress 
quantum ones. At very high temperature the classical 
term is dominant. It is clear that Ea. l|3.6[l is an expan- 
sion in powers of f3J/s — Js/T, so it is valid as far as 
T> Js. 

To calculate the partition function we use 
Ecis. lf2~(3|) . lf2~7|) where H comes from Eg . (|3~j)|) . Inte- 
gration over spherical angles of ni and n 2 leads to the 



\nZ = ln[(2ws + l)(2s + 1)] + In 
B(l + w)£ 1 



sinh(£w) 



2s 12s 2 



2 -(TS^ 2 -3^ 2 + 15^ 
-2w 2 £ 2 + &Buji - 2£ 2 ± — + 3(1 + Lofi 2 



LO 

9(1 +co) 2 3B 2 ^(1+uj) 2 



(3.7) 



where £ = (3 J and B = coth(w£) — l/(u;£) is the 
Langevin function. We can obtain the internal energy 
(U = -dhxZ/d/3) and heat capacity (C = dU/dT) in 
this approximation. Quasiclassical expansions for the in- 
ternal energy and heat capacity of this model are shown 
in Figs. pi2l3l4| ) for {S u s 2 )=(5/2, 2). They show how 
different contribution of classical, first and second order 
corrections sum up to the reasonable result. We have 
also plotted the exact solution of ferrimagnetic gas model 
for comparison. The corresponding values for ferrimag- 
netic chain of interacting molecules have been also plot- 
ted which will be discussed in next sections. In all of 
these plots the horizontal axis is scaled to dimensionless 
parameter t — T / Js 2 — T / J. Thus for large spins, it cov- 
ers even the intermediate and low temperature regions. 
In the present case the convergence range of cumulant 
expansion reads to t > 1/s = 1/2. 

A good agreement is observed between second order 
(0(l/s 2 )) cumulant expansion and exact result for inter- 
nal energy of F coupling which is shown in Fig.Q. For 
t < 1/2 it shows large deviation from exact result and di- 
verges as t — ► 0. The reason for divergence is as follows. 
The internal energy up to second order approximation 
(0(l/s 2 )) would be, 



U 
Is 2 



\ujB + ^^{Ai + B) 



2s 



1 



[( =F 2Blo - 6lo + 15A - Alo 2 + 6Au - 4 
GA _ „ „(l+w) 2 



12s 2 

+ — +6(1 + cj) 2 -9 



A- 3(1 +lo) 2 B 2 )S, 
6B 



=F Auj - 3(1 + uj) 2 BA)^ 2 + 15B + 6Blo 



(3.8) 



where A = dB/d£. As T — > 0, U diverges since £ — > 
oo. The coefficient of linear term in £ comes from both 
Tt^ and Ti^ where they cancel each other as £ — > oo. 
Whereas the quadratic term £ 2 needs a contribution from 
Ti^ to be regular at £ — > oo, which is not present in 
second order approximation. The heat capacity for F 
coupling of ferrimagnetic molecule is shown in Fig- (O- 
The classical and first order expansion are far from the 
exact result, however the correction of second order term 
makes an agreement for t > 1/2. 
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FIG. 1: Internal energy per molecule of the ferrimagnetic 
(Si — 5/2, S2 = 2) gas model with ferromagnetic (F) coupling. 
Exact result of gas model (solid line) and different orders of 
cumulant expansion for gas model and chain of interacting 
molecules. 



FIG. 2: The heat capacity per molecule of the ferrimagnetic 
(Si = 5/2, S2 = 2) gas model with ferromagnetic (F) coupling. 
Exact result of gas model (solid line) and different orders of 
cumulant expansion for gas model and chain of interacting 
molecules. 



IV. FERRIMAGNETIC SPIN CHAIN 



For AF coupling of gas model the normalized inter- 
nal energy (U / J) is plotted in Fig. (3). It is surprising 
that the second order expansion fits very well on the 
exact result even for very small values of i. A simi- 
lar behavior can be seen in Fig. 10} for heat capacity. 
To understand the different behavior of cumulant ex- 
pansion for F and AF coupling, let us come back to 
the energy levels. For F coupling the ground state has 
S]P = (Si + s 2 ) and any excited state (E n ) can be ob- 
tained by S[ n ^ = (Si + s 2 ) — n, n = 1, . . . , 2s 2 , (Si > s 2 ). 
In the spectrum of F case the difference of energy levels 
(AE^ = Ei — Ei-i = JS^ 1 ') decreases by increasing 
Ei . So it is like a quantum system. Where the converse is 
true for AF case. Moreover the absolute value of AE for 
the lowest states of F coupling is larger than AF case, i.e. 
AE^(F) - AE^(AF) = J(2s 2 - 1). This means that 
cumulant expansion works better for AF case. In this 
respect we should note that a cumulant expansion sim- 
ilar to a high temperature series expansion starts from 
T — > oo, where the probability of all states are equal. 
Then different order of expansion are responsible to re- 
cover the non-equal probability of states. This is cru- 
cial when the low energy spectrum is not dense (like F 
case) where a big difference exists between the occupation 
probability of lowest levels. Thus we expect to observe 
stronger quantum effects in the ferromagnetic case. 



We will now consider an interacting system of ferri- 
magnetic molecules. As it was noted in the introduction 
such systems can be synthesized to be considered effec- 
tively as one dimensional models (chains)^, since the in- 
terchain coupling is very small compared to intra-chain 
ones. In this case the Hamiltonian of a ferrimagnetic 
chain of molecules is: 

JV JV 
2 ' 2 

H = TJ2 J 2i-l,2jS 2 i-l ' S 2j (4.1) 

i,j 

It is supposed that larger spins (S) sit on the odd- 
numbered sites and smaller ones (s) on the even sites. 
The exchange coupling (J2j-i,2j) exists between any dif- 
ferent spins. Our approach is general to cover long-range 
cases, however, later we will consider nearest neighbor 
coupling to obtain physical quantities. 

To find an expansion for the physical quantities we 
should first calculate the semiclassical effective Hamil- 
tonian defined in Eci. (|2.9[l . Again, we will consider the 
cumulant expansion up to second order corrections of 1/s 
which comes from TL ^ . Let us first calculate the classical 
term which is the simplest one. 

(H) c = (H) =T^ti, 2] (S„ • s 2i ) 

= Tw^J 2! -i,2jn2 i -i-n2 3 (4.2) 

i,3 
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FIG. 3: Internal energy per molecule of the ferrimagnetic 
(Si — 5/2, S2 = 2) gas model with antiferromagnetic (AF) 
coupling. Exact result of gas model (solid line) and differ- 
ent orders of cumulant expansion for gas model and chain of 
interacting molecules. 



FIG. 4: The heat capacity per molecule of the ferrimagnetic 
(Si = 5/2, s-2 — 2) gas model with antiferromagnetic (AF) 
coupling. Exact result of gas model (solid line) and differ- 
ent orders of cumulant expansion for gas model and chain of 
interacting molecules. 



where as before J 2 i-i )2 j — J2i-i,2jS and Ea. l|3.3[l has 
been used. The first correction (0(1/ s)) comes from the 
following term. 

(HH) C = ^2 ^2i-l,2j^2/-l,2fe((S2i-l • S2j)(S2;_l • S 2 fc)) C 
i,j,l,k 

(4.3) 

There are different sequences of i,j,l and k in ((S2i-i ■ 
s 2j)(S2;-i ■ 82k)) where indices can be in the same site, 
neighboring sites or apart from each other. We have sum- 
marized different cases in appendix which finally leads to 
the following result for TiS 1 ^. 

= ^E^- 1 ,2 J (l-n2 l - 1 -n 2j ) 2 

(n 2j ■ n 2 / - (n 2 i_i • n 2j -)(n 2i _i • n 2 /)) 
fa? v- j J 

2-1 ^i-l^j J U~l.2j X 

(n 2 i_i • n 2i _i - (n 2l _i • n 2i )(n 2 ;_i • n 2j )) 

(4.4) 

As it is obvious from Ea. (|4.4|l . the first correction which 
contains quantum effects is not of Hcisenberg type. 
Moreover the second and third term in Eq. I|4.4|l contains 
coupled interaction of 3 sites, i.e. (n 2 j_i ■n 2 j)(n 2 j_i -n2z) 



and (n 2 ,_i ■ n 2 j)(n 2 ;_i ■ n 2 j). This is different from two 
sites interaction of the original Hamiltonian, Ea. H4.1|) . It 
is the price of working with a classical Hamiltonian in- 
stead of the original quantum form. The last two terms 
of Ea. l|4.4fl give the information of 3 sites correlation. 
We will come back to this point later when comparing 
the results of gas model with a chain of interacting fer- 
rimagnetic molecules. The third term of cumulant ex- 
pansion is (HHH) C , where the cumulant of 6 operators 
should be calculated. Each cumulant is composed of such 
terms, ((S 2i _i • s 2j -)(S 2 /-i • s 2fc )(S 2m _i • s 2 „)) c - We have 
shown in appendix the different possible sequences of 
(i,j,l,k,m,n) which give nonzero cumulant. Summing 
up all nonzero cumulants, the second order semiclassical 
correction is 



uj(3 2 

i : j,l : k 

- ^ V J J J e 

12s 2 ' ' J 2i-l,2j J 2k-l,2j J 2l-l,2j^3 
T / ](jji-l,2i J2i-l,2iQ4 + W^2i-l, 2 i^2i-l,2j@5) 

=F ^iE^-i^ee (4.5) 
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where 0; 



. , 6 are defined below. 



V. RESULTS AND DISCUSIONS 



©l = (n 2 j • n 2fc _i) - (n 2l _i • n 2j )(n 2j _i • n 2fc _i) 
-(n 2 j • n 2/ )(n 2fe _ 1 • n 2/ ) 
+(n 2 i_i ■ n 2j -)(n 2i _i ■ n 2/ )(n 2fc _ 1 ■ n 2/ ) 

+ ^((n 2i -i • n 2 fe_i)(n 2 j ■ n 2; ) 
-(n 2i _i • n 2 ;)(n 2j - • n 2 fe_i)), 
6 2 = -4(n+_ 1 • n 2j )(n 2i _! • n 2 0(nii_i • n 2fe ) 
= (n 2i _i • n 2i )(n 2i _i ■ n 2 ;)(n 2l _i • n 2fc ) 

-^[(n 2i _i • n 2j )(n 2/ • n 2fc ) 

+ (n 2l _i • n 2; )(n 2j - • n 2fc ) 
+ (n 2l _i • n 2fc )(n 2 j • n 2/ )], 

9 3 = -4(n 2j _i ■ n^.)(n2fc_i • n 2j )(n 2 ;_i ■ n^-) 
= (n 2i _i • n 2j -)(n 2fe _i • n 2j )(n 2/ _i • n 2j ) 

~[(n.2i-i ■ n 2j -)(n 2fe _i • n 2; _i) 

+(n 2fe _i ■ n 2j )(n 2i _i ■ n 2 ;_ x ) 
+ (n 2 ;-i ■ n2j)(n 2i _i -n^i)], 

9 4 = -{(n 2j _i ■ n 2; ) (1 - (n 2i _i ■ n 2i )) 2 

+2[(n 2 j • n 2/ ) + (n 2i _i • n 2i ) 2 (n 2 . t _i ■ n 2/ ) 
-(n 2j _! ■ n 2j )(n 2j • n 2/ ) 
-(n 2i _! ■ n 2i )(n 2i _i ■ n 2/ )]}, 

9 5 = -{(n 2; _i • n 2j -) (1 - (n 2l _i • n 2j )) 2 

+2[(n 2i _i • n 2 /_i) - (n 2i _i • n 2j -)(n 2i _i • n 2 /_i) 
-(n 2i _i • n 2j -)(n 2 ;_i • n 2j ) 
+ (n 2; _i • n 2j -)(n 2 j_i • n 2j ) 2 ]}, 

9 6 = (n 2i _i • n 2j -) (1 - (n 2l _i • n 2j )) 2 . (4.6) 

Let us draw your attention to two properties of Eq. 14.5J1 . 

If we factorize J — Js 2 from each term they are in (^f) 2 
order. However there are extra prefactors of powers of 
(1/s) for some terms which make them less important for 
large s. This is a general property which is also valid for 
other orders of cumulant expansion. Generally, in the cu- 
mulant expansion the n-th order term, Tv- n \ can be con- 
sidered of J(^f) n order. The second point comes again 
from non-Hcisenberg type of interaction in O j . There are 
some terms which contain coupled interaction of 4 differ- 
ent sites, like (n 2i _i • n 2j )(n 2i _i • n 2; )(n 2fe _i • n 2/ ) in 9i. 
Such terms give the information of 4-point correlation 
function. Thus the n-th order semiclassical Hamiltonian, 
, recovers the information of (n + 2)-sites correlation. 
However we found that for nearest neighbor interaction 
(Jij = JSi.j±i) the expectation value of contains 
only nearest neighbor interaction. Thus, it contains the 
correlation of two sites the same as TiS 1 ^ for gas model. 
This means the first order approximation of ferrimagnetic 
gas model (H^ +TiP^) would give the same approxima- 
tion of first order one dimensional interacting molecules. 
This will be shown in next section. 



In this section we will present some physical quantities 
of a ferrimagnetic chain. Our results of previous section 
are general for any long-rang interaction of Jij. How- 
ever we consider nearest neighbor case, Jy = J<£j,j±i, to 
obtain internal energy and heat capacity. Because most 
of synthesized materials behave as nearest neighbor in- 
teracting molecules, moreover a comparison with other 
results is in this case. 

Similar to the case of gas model the implementation 
of Eas. (|2.6|) and (|2.7[1 where TC is approximated by 

H (o) + H {i) + n m leads t0 



z = z [i -p(H {1> +n {2) ) 



2! 



(!)1 2 \ 



and 



4?r 
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V 4tt 
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To 2nd order approximation the partition function is 

n ( sinh(u£) x A : 



Z 



(2ws + l)~(2s + l) 



B^ t (l+u)B£ 



2s 



+N 



<f / 2Bw 6B 2uj\ 

4?v ~W + ~e~T) 



+ 



6s 2 



1 3B B 2 



?2,.,2 



+ ]y | £ n , ,.\2 fn '''' a- ,,2 

8s 2 



(1+ujy 2 - 1- NB - 

1 ^UJ ' 



cues 2 



4s 2 



(5.1) 



where we have used some integrals of Ref . . Taking 
derivative of In Z to (3 gives the internal energy 



U 
J 



luB 



1+UJ 

2s 



(M + B) 



+ — [(2AB£ 2 + 2B 2 £) (-2 - 2lo 2 - 'ioj) 

+ (A( + 5)(15 + 6co + - - -(1 + c^) 2 ) 
+2£_ (-3u - 2 - 2uj 2 + 3(1 + uj) 2 ) 

^u{Ae + 2BQ] 



(5.2) 



That's simple to calculate the expression of heat capacity, 
C = dU/dT, which is not presented here. 

We have plotted in Figs. (^13) the internal energy per 
molecule (2U/NJ) of (5i = 5/2, s 2 = 2) ferrimagnetic 
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chain ( line) for F and AF coupling respectively. 

As far as we know there is no other result for large spin 
systems to compare with. However we compare them 
with the corresponding values of different orders of gas 
model. We observe that the zero (classical) and first 
order expansion is the same for gas model and interacting 
chain. Because the 1st order expansion (H^ + HS 1 ^ 1 ) 
contains information of two sites correlation. This can 
be shown easily by the following equation 

R J2 N 

= ^-£((l- ni -n i+1 ) 2 ) 

l — l 

BP N 
_^_ w(w + l)^ (1 _ (n .. n . +l) 2 ) . 

i=i 

(5.3) 

As far as the correlation of two sites is concerned the 
behavior of a single molecule and a chain of molecules 
with nearest neighbor interaction is the same. The in- 
set of different behavior for chain and non-interacting 
molecules comes from the correlation of 3 sites. Such 
terms exist in Ti.^ and contribute in 1/s 2 correction. We 
see the difference by 1/s 2 expansion in Fig.Q). This de- 
viation is more clear for low temperature regime where 
quantum effects are important. But for temperatures 
greater than classical energy scale, T > JS1S2, there is 
a good agreement between gas model and ferrimagnetic 
chain. This shows that the correlation of more than two 
sites is important for low temperatures, t < u). In other 
words, at moderate and higher temperature an ideal gas 
of molecules represents a chain of ferrimagnets very well. 

We have also plotted in Figs.© and (0J the heat ca- 
pacity of ferrimagnetic chain for F and AF coupling re- 
spectively. Similar to the behavior of internal energy, the 
classical and 1st order (1/s) expansion gives the same 
results of gas model. The 2nd order correction (1/s 2 ) 
makes the difference between chain and single molecule. 
However a single molecule represents well the behavior 
of a chain. This justifies long range correlations can be 
neglected for moderate and high temperature regime of 
ferrimagnetic chain. This is in agreement with previ- 
ous conclusion that we can replace the results of the gas 
model for a chain of interacting ferrimagnetic molecules 
where temperature is greater than classical energy scale, 

t > OJ. 

Heat capacity of both Figs.J21 decreases for high t 
which is the sign of antiferromagnetic behavior in the up- 
per part of spectrum of a ferrimagnetic chain. We also see 
a Schottky-like peak in C which is the result of ferromag- 
netic to antiferromagnetic crossover. However this peak 
is in the region where cumulant expansion is not neces- 
sarily convergent. Then our plots for t < 1 might not be 
reliable, although we know from other arguments^ this 
peak exists. 

To see the effect of spin inhomogeneity we have also 
considered the case of S\ = 4, S2 = 2 (oj = 2) for both 
gas model and ferrimagnetic chain which have not plot- 



ted here. As mentioned before the first order correction 
contains information of two sites correlation which makes 
no difference between gas model and a chain of interact- 
ing molecules. We have found that quantum effects are 
more pronounced for larger w in F coupling both in the 
internal energy and heat capacity. The second order cu- 
mulant expansion of heat capacity for larger oj reaches 
the exact value for temperature higher than the case of 
Fig. • Although we expect a better agreement for lower 
temperature by increasing spin magnitude, this is not the 
case if we increase one of them. Specially for F coupling, 
increasing oj makes larger energy difference for lower lev- 
els which makes the situation worse. 



1.2 




t=T/Js 2 

FIG. 5: Heat capacity per molecule of ferrimagnetic chain of 
(3/2, 1) with AF coupling. Second order cumulant expansion 
(CE), Monte-Carlo (MC) simulation and spin wave theory 
(SWT). 

For the AF coupling we have observed almost no dif- 
ferent behavior comparing Figs. (J3I ^Ql (w = 5/4) with the 
results of oj = 2. In both groups, the second order cumu- 
lant expansion reaches the exact result of gas model at 
very low temperature. However this will not remain the 
same for large oj. We have increased oj to S1/S2 = 7/2 
and observed good agreement for a little higher temper- 
ature. But still fairly well results of cumulants in the 
convergent region t > 1/2. 

To the best of our knowledge there is no high tem- 
perature expansion for the thermodynamic properties of 
ferrimagnets in a closed form in the literature. The large 
S spin wave approximations^*^ are imposed to an extra 
constraint to be valid for moderate temperature. Because 
it is valid for low temperature (by definition) where only 
the low energy spectrum has the dominant effect. Adding 
an extra constraint needs a self consistent numerical solu- 
tion to obtain thermodynamic properties. In comparison, 
our results are analytic (in a closed form) and general 



8 



in the sense that it can be applied to any ferrimagnetic 
system for arbitrary exchange couplings. The effective 
Hamiltonian derived here are in terms of arbitrary Jjj 
which covers long range interactions as well as any lattice 
structures. 

However, we finally present the result of second or- 
der cumulant expansion (CE), spin wave theory (SWT) 
and Monte-Carlo (MC) simulation^ for heat capacity of 
(S\ = 3/2, s 2 — 1) ferrimagnetic chain in Fig.JSJ. Al- 
though cumulant expansion is supposed to work well for 
large spin systems we present our results for (3/2, 1) fer- 
rimagnets to compare with available results. The qual- 
itative behavior of CE is the same as MC simulation 
but large deviation is observed for moderate tempera- 
ture. The main reason for this deviation is the strong 
quantum nature of small spins. However, much better 
agreement is expected for larger spins. 

As an out look, this work can be generalized to ladder 
geometry. It is known that the zero temperature behavior 
of ferrimagnetic ladders^! is different from homogeneous 
spin counterparts 2 ^. It is interesting to see the differ- 
ence for finite temperature, for instance the evolution of 
magnetization plateau x 23 i 24 . 
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APPENDIX A: THE CUMULANT OF FOUR AND 
SIX OPERATORS FOR A SPIN CHAIN 

In order to achieve the first order cumulant i.e. ((S2«-i ■ 
S2j)(S2i— l • S2fe)) c , there are three nonzero sequences of 
i, j, / and k as follows 

• Nonzero cumulant for i = I and j = k. These are 
listed in table |U 

• For i = I and j =/= fc, the cumulant is 
((S2i-i • S2j)(S2/-i • S2j)) c and its value equals to 
2cj 2 s 3 (n 2i -i • n 2 F J -)(n 2 ;-i • n^) 

• In the case of j = k and j ^ I the cumulant 
is ((S 2l -i • s 2 j)(S 2 i-i • s 2fe )) c which is equal to 
2ws 3 (n+_ 1 • n 2j )(n 2i _ 1 • n 2; ) 

In order to achieve the second order cumulant of the 
spin chain i.e. ((S 2i _i • s 2i )(S 2i _i • s 2fe )(S 2m _i • s 2 „)) c , 
there are four nonzero sequences of i, j, I, k, m and n as 
follows 

1. Three operators sit at a site and the others at an- 
other site, the cumulant is ((S 2 j_i • s 2 j)(S 2j _ 1 • 
s 2j)(S 2 i_i • s 2 j)) c . There are only 15 different se- 
quences which are nonzero. 



TABLE II: Nonzero cumulant of six operators in a chain of 
interacting molecules (Ea. J4.ip - ). The first column shows the 
position of different spins. For instance 135, 24, 6 shows that 
Xi , Xi and X5 sit at a single site, X2 and £4 at the other site 
and xe sits at a separate one. (A) c = ((XiX2)(X3X4,)(XsXe)) c - 









X2 


x 3 




X4 


x 5 




xe 


(A) c 


135,24,6 


"->2i- 


1 


_|_ 


S21- 


1 


S 2j 


&2i- 


1 


s 2l 






Sai- 


1 


4 


&2i- 


1 


S 2j 


&2i- 


1 


a 21 


9 A 




Sti- 


1 


i j 


S%i- 


1 


S 2j 


S21- 


1 


S21 


Q 




Sti- 


1 




S%i- 


1 


Sfj 


S21- 


1 


s 2l 


-2lus 4 


135,26,4 


&2i- 


1 


4 


S%i- 


1 


821 


Sii- 


1 


S 2j 


lo 2 s 4 




&2i- 


1 


4 


&2i- 


1 


a 21 


&2i- 


1 


S 2j 


9 A 




Sti- 


1 


4 


sii- 


1 


S21 


&2i- 


1 


S 2j 


-4ws 3 




Sti- 


1 


• s 2j 


S21- 


1 


a 21 


"2i- 


1 


S 2j 


-2us 4 


135,46,2 




1 


six 


S%,i- 


1 


4 


&2i- 


1 


S 2j 


4lj 2 s 4 




Sti- 


1 


all 


&2i- 


1 


4 


Sii- 


1 


S 2j 






Sti- 


1 


s z 2i 


S%i- 


1 


4 


"2i- 


1 


S 2j 






si- 


1 


S Z 2l 


S%i- 


1 




&2i- 


1 


S 2j 


-2lus a 


246, 13, 5 


c+ 

b 2i- 


1 




&2i- 


1 


4 


S2k- 


-1 


S 2j 


A 2 \ 

A10 s 




0+ 
b 2i- 


1 


: 


$2i- 


1 




S^k- 


-1 


S 2j 


a 2 4 

4u> s 




sii- 


1 


4- 

4 


&2i- 


1 


S 2j 


•S'lfc- 


-1 


S 2j 


a 2 3 

-4u> s 




LJ 2i- 


1 


at 
"2j 


^2i- 


1 


So- 


C- 

> - J 2k- 


-1 


"2 3 


-2lo 3 s 4 


246, 15,3 


Sti- 


1 


s+ 

S 2j 


qZ 

°2k- 


-1 


S 2j 


&2i- 


1 


S 2j 


4c^V 




$2i- 


1 


a 2 j 


Sik- 


-1 


• S 2j 


&2i- 


1 


S 2j 


4^ 2 s 4 




si- 


1 


S 2j 


qZ 

°2k- 


-1 


o z 
S 2j 


"2i- 


1 


S 2j 


-4lj 2 s 3 




&2i- 


1 


S+ 


S2k- 


-1 


„z 
S 2j 


sii- 


1 


S 2j 


-2uj a s 4 


246, 35, 1 


°2k- 


-1 


S 2j 


D 2i- 


1 


o z 
S 2j 


"2i- 


1 


S 2j 


4^ 2 s 4 




°2k- 


-1 


S 2j 


&2i- 


1 


S 2j 


^21- 


1 


S 2j 


4^V 




qZ 

°2k- 


-1 


S 2j 


°2i- 


1 


Q Z 

«2j 


&2i- 


1 


S 2j 


-4^ 2 s 3 




S2k- 


-1 


S+ 


S2i- 


1 


S 2j 


Sii- 


1 


S 2J 


-2^ 3 S 4 



2. Three operators sit at a site, two operators sit at 
another one and the sixth one sits at a separate 
site, the nonzero cumulants are listed in table [H] 
At the first column of the table |nj we can observe 
operators which can be placed at the same site, (i.e 
at the first four row of the table i = I — k, j = I 
and k 7^ n) 

3. When two operators sit at a site, the another two 
at another site and the fifth and sixth one sit at 
two different sites, the cumulant results to 

4w 2 s 4 (2(n+_ 1 • n+ )(n 2 " i _ 1 • n 2i )(n 2i; _i • n^-) 

+ (nj-i ' n 2/ )(n^-i ' n2j)( n 2fc-i ■ ™2j) 
+ 2 ( n i-i ■ n2/)(n 2fe _i ■ n+ )(n 2 " i _ 1 ■ n^) 

+ ( n 2fc-l ' ^tj)( n 2i-l ■ n 2j)( n 2i-l ' n 2/))- 

4. Finally, if three operators sit at a site and other 
ones sit at three different sites, then there are two 
nonzero terms as follows 
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(a) If i = I = m, j, k and n are not equal, the 
value of cumulant is 

-2ws 4 (n+_ 1 • n 2i )(n 2i _i • n 2 /)(n 2i _ 1 ■ n 2fe ). 

(b) If j = k = n, i, I and m are not equal, the 
value of cumulant is 



2w 3 s 4 (n 2j _i • n+ )(n 2 fe_i • n 2i )(n 2; _i • n 



2j> 



Note that, the second order cumulant expansion is the 
sum of above mentioned cases. 
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